Polarizabilities of Germanium Clusters 
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Polarizabilities of Ge„ clusters with 2 to 25 atoms are 
calculated using coupled-perturbation Hartree-Fock (CPHF) 
and finite field (FF) method within density functional the- 
ory. The polarizabilities of the Ge„ clusters increase rapidly 
in the size range of 2 to 5 atoms and then fluctuate around 
the bulk value. The polarizabilities are sensitively dependent 
on the cluster geometries and electronic structures. The large 
HOMO-LUMO gap may lead to the small polarizability. As 
compared with the compact structure and diamond structure, 
the prolate cluster structure corresponds to a larger polariz- 
ability. 
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In the past two decades, the structural and electronic 
properties of semiconductor clusters have been exten- 
sively studied because of their fundamental, interest and 
potential application in nanoelectronicsoc3. The small 
semiconductor clusters are well understood up to 10 
atoms. But our knowledge for larger clusters is still quite 
limited. The polarizability is one of the most impor- 
tant quantities of the clusters, which can yield the static 
dielectric constant in the bulk limit through Classius- 
Mosotti relation. On the other hand, the polarizability 
can provide some information on the bonding and geo- 
metrical features of the clusters. Thus, comprehensively 
understanding of the polarizabilities from theoretical cal- 
culations— is important in cluster science. The jellium 
modeHEZl was successfully applied to study the polariz- 
ability of large metallic clusters. But for semiconductor 
clusters, the bonding and geometrical effects are not in- 
corporated in the jellium model. Alternatively, the ab 
initio calculations based on quantum chemistry meth- 
ods are needed. There were only few previous ab fflitip 
calculations on the polarizabilities of the clustersETE 3 ! 
Especially, for the germanium clusters, there is only one 
attempt to calculate,the polarizabilities and the cluster 
size is quite limitednS. _ 

In our previous studiescj, the geometries of Ge n clus- 
ters (n = 2 — 25) have been obtained by density func- 
tional DMalj-calculationsa incorporated with a genetic 
algorithmEj El It was found that the Ge„ clusters follow 
a prolate growth pattern with n > 13. The stacked lay- 
ered structure and the compact structure compete with 
each other in intermediate size range. Based on the previ- 
ously optimized low-energy structures, in this paper, we 



will study the polarizabilities of these clusters using the 
analytically coupled-perturbation Hartree-Fock (CPHF) 
and numerically finite field (FF) methods within density 
functional theory. We aim to explore the size dependence 
of the polarizability, the influence of the atomic and elec- 
tronic structures on the polarizabilities of the Ge„ clus- 
ter. 

It is well known that electron correlation plays a pri- 
mary role in determining molecular polarizabilities. DFT 
treatment has been proven to make significant improve- 
ment to Hartree-Fock results in moleculax-.polarizabili- 
ties and hyperpolarizabilities calculationso'EJ. Within 
DFT framework, B3LYP functional considers the hy- 
brid betaseen the Hartree-Fock exchange-and Kohn-Sham 
orbitalsEJO, while LANL2DZ basis setc3 can give a good 
description of the bonding and geometrical features of 
heavy atoms. Thus, B3LYP /LANL2DZ scheme is ex- 
pected to well describe the cluster polarizabilities at ac- 
ceptable computational cost. Here, all these calculations 
are performed alB3LYP/LANL2DZ level by using Gaus- 
sian ^ packaged. 

The energy in an external electric field can be ex- 
panded as 

1 

— ( 

2 



e{f) = e{o) - - /•',/••, - 



(1) 



where E(0) is the energy without the external field, Fi 
are components of the applied field. The dipole moment 
Hi and the polarizability ay are defined as: 
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Within Kurtz's finite field methodEl, the dipole moment 
and polarizability in a uniform field, can be derived from 
the follow equations. 

IHFi = ~\ [E(F) - E(—Fi)] - i [E(2F l ) - B(-2F,)] 



a u F t = -E(0) - - [E(Fi) + E{~F i: 
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[E(2Fi) - E(-2Fi)] 
(5) 
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To calculate of the dipole n and the polarizability a, 
at least 13 self-consistent field (SCF) runs are neces- 
sary with the field strengths ±Fi and ±2Fi (i = x, y, z). 
One of the most severe problems in finite-ficM method 
is the choice of an appropriate field strengtheH. Sim et 
al.C3 have assessed the numerical accuracy against dif- 
ferent field values and concluded that stable linear and 
nonlinear polarizabilities can be obtained when F equals 
0.001 ~ 0.005 a.u. In this work, the external field is 
added along x, y, z with the magnitude 0.001 a.u and a 
tighter SCF convergence criterion of 10 -9 is adopted. 

The measured data in experiments are usually the av- 
erage polarizabilities, which can be obtained by 

< a >= -(a xx + a yy + a zz ) (6) 

The optimized structures for Gen_25 clusters at 
B3LYP /LANL2DZ level are shown in Fig.l. Similar 
equilibrium configurations are found by DMol and Gaus- 
sian98. Therefore, the details for the structures will not 
be discussed here. 
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FIG. 1. Lowest energy structures for Ge n (n = 11 — 25) 
clusters. 



Table I. The dipole moments /x(Debyes) , polarizabilities 
a(A2iatom) of— Si2-io clusters :a a : present results, a b : 
Ref.E3,a c : RefB, a d : RefE3, (in bracket) :Ref.E3. 
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0(0) 


4.97 


6.29 
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0.32(0.33) 


4.56 


5.22 
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0(0) 


4.48 


5.07 
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0(0) 


4.87 


4.81 
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0.01(0.0) 


4.62 


4.46 
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0(0) 


4.60 


4.37 
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0.23(0.0) 


4.83 


4.52 






9 


0.17(0.36) 


4.73 


4.38 


4.46 


3.0 


10 


0.29(0.69) 


4.55 


4.31 


4.65 


5.50 



To check the validity of current method, we firstly cal- 
culate the dipole polarizability of small silicon clusters. 
Table I compares our calculations with previous theoret- 
ical and experimental results. The current theoretical 
dipole moments fox. the clusters Si2-7 are in agreement 
with Vasiliev et ac3. But our results for the dipole mo- 
ment of Sis-it) clusters are different due to the substantial 
difference in geometrical configurations. The present geo- 
metric structures for Si clusters are consistent pHiith those 
obtained by Shvartsburg et al. and Li et al£3t3. Our 
calculated polarizabilities seem to be overestimated as 
compared with Vasiliev et ae3, but consistent with Deng 
et aE3. Moreover, the present polarizabilities are sensi- 
tively dependent on the cluster size and oscillate with the 
cluster size. While the polarizabilities in Vasiliev's work 
tend to decrease with the increase of the size. Exper- 
iments show that the polarizabilities fluctuate with the 
cluster size. 

Table II. The dipole moments /i , polarizabilities a (A 3 / 
atom), binding energies Ef, (eV), HOMO-LUMO gaps A 
(eV) of Ge„ clusters: ^cPHF(Debyes): XjPHF method; 
/ipp (a.u.): FF method; in bracket: Ref.c3. 
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MCPHF (PRL) 


Mff 
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E b 
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0(0) 





4.10(6.67) 


1.69 


0.31 
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0.61(0.43) 


0.241 


5.07(5.89) 


2.84 


1.32 
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0(0) 





5.14(5.45) 


3.23 


1.14 
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0(0) 





5.52(5.15) 


3.32 


1.48 
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0.15(0.14) 





5.36(4.88) 


3.41 


1.39 
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0(0) 





5.27(4.70) 


3.50 


1.36 
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0.55(0) 


0.22 


5.47(4.99) 


3.44 


1.10 
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0.12(0.28) 


0.05 


5.39(4.74) 


3.48 


1.22 


10 


0.56(0.68) 


0.22 


5.17(4.66) 


3.59 


1.31 


11 


1.35 


0.53 


5.24 


3.53 


0.99 


12 


1.58 


0.62 


5.33 


3.49 


1.05 


13 


0.82 


0.32 


5.45 


3.57 


0.98 


14 


1.30 


0.51 


5.41 


3.61 


1.20 


15 


0.18 


0.07 


5.44 


3.57 


0.80 


16 


0.74 


0.29 


5.39 


3.61 


1.17 


17 


1.11 


0.24 


5.52 


3.59 


0.94 


18 


0.61 


2.75 


5.61 


3.58 


0.90 


19 


0.46 


0.18 


5.41 


3.63 


0.70 


20 


1.06 


0.46 


5.51 


3.61 


0.89 


21 


1.70 


0.67 


5.54 


3.62 


0.94 


22 


0.51 


0.20 


5.45 


3.63 


0.81 


23 


2.89 


1.14 


5.62 


3.62 


0.82 


24 


2.52 


0.99 


5.69 


3.64 


0.62 


25 


0.52 


0.21 


5.63 


3.61 


0.70 
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Table II gives the dipole moments, average polarizabil- 
ities, binding energies and HOMO-LUMO gap as func- 
tions of the cluster size. The dipole moments calculated 
with CPHF method and FF method are compared and 
agree well with each other. The dipole moment reflects 
the symmetry of the geometrical structure: the smaller 
dipole moment corresponds to the higher symmetry. The 
dipole moments for the clusters with n = 2,4,5,7 are 
nearly zero, corresponding to the high symmetry in these 
clusters. For Ge4, the lowest energy structure is T>2h 
rhombus, while they arc trigonal bipyramid with D 3 ^ and 
pentagonal bipyramid (D§h) for Ges and Gey. In the case 
of Geg , Geg and Geio, capped pentagonal bipyramid, 
tricapped trigonal prism and tetracapped trigonal prism 
have favorable energy, respectively. The dipole moments 
and the polarizabilities of small Ge clusters are compared 
with Vasiliev et alEB and our calculations seem to be a 
little overestimated. 



o 




& 4.5 



> 
< 



10 15 

Cluster size 

FIG. 2. The polarizabilities of Ge clusters are plotted as a 
function of the cluster size for Ge n . 

Fig. 2 gives the polarizabilities of Ge„ clusters as a func- 
tion of cluster size. The polarizabilities of small clusters 
increase rapidly with cluster size in the size range of 2 
to 5 atoms. Then, the polarizabilities fluctuate around 
a = 5.4A 3 /atom. In the light of Clausius-Mossotti rela- 
tion, 



eM + 2 



(7) 



the polarizability of the bulk Ge is 4.5A 3 /atom. The 
present average polarizability of the clusters is larger than 
that of the bulk. Previous studies have shown tha^-jthe 
prolate structure contributes a large polarizabilitytBEj. 
The prolate configuration is preferred for medium-sized 
Ge clusters, which occupies larger distortion and lower 
symmetry in comparison with diamond structure. Thus, 
the polarizability of the clusters is larger than that of the 
bulk. It is worthy to note, that the present results is quite 
different from Vasilier'sES. Compared with their method, 
the present scheme incorporate the electron correlation 
in the calculation of the polarizability of small clusters. 

We further discuss the relationship between the polar- 



izability and the electronic structures of clusters. Fig. 3 
plots the polarizabilities as a function of the HOMO- 
LUMO (highest occupied molecular orbital (HOMO) and 
lowest occupied molecular orbital (LUMO)) energy gaps. 
As shown in Fig. 3, the large polarizability generally cor- 
responds to the small HOMO-LUMO gap. The clusters 
with n = 8, 11, 13, 15 have relatively smaller polarizabil- 
ities in comparison with their neighboring size. For ex- 
ample, the HOMO-LUMO gap for Ge 8 is l.lOeV and the 
polarizability is just 5.17 A 3 /atom. The polarizabilities 
decrease with the increase of the HOMO-LUMO gap ex- 
cept for the clusters Ge2, Ges and Geig_22- p^rUf can be 
easily rationalized using the two-level modelEirE3, 



a 



(8) 



where is the transition dipole moment from the 
ground state to the first dipole-allowed excited state and 
At the corresponding transition energy. Approximately, 
At can be replaced with HOMO-LUMO energy gap A. 
From this model, a increases with decreasing A , con- 
sistent with our calculated trend for most Ge clusters. 
However, a is not a simply inverse proportion relation 
to A, since the two quantities are dependent each other. 
Small A tends to generate large fj,f. In addition, $ de- 
pends on some other characters, such as selection rule. 
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FIG. 3. Polarizabilities of Ge clusters are plotted against 
the HOMO-LUMO gap. 

However, the polarizability of the cluster is not only de- 
pendent on HOMO-LUMO gap but also closely related 
to geometrical characteristics. For example, the HOMO- 
LUMO gap of Geis is larger than that of Gei9, while the 
polarizability of the former is larger than that of the lat- 
ter. This can be attributed to their structural difference. 
For Geis, the optimized structure is two interpenetrate 
pentagonal connected with a bicapped square antiprism 
Geio subunit, while a more compact cage-like geometry 
is favorable to Geig. The compact structures have rela- 
tively less and shorter bonds, which leads to the valence 
electrons binding tighter. Thus, a smaller volume is ob- 
tained in the compact structure, which causes a smaller 
polarizability for Geig. Similarly, the different polariz- 
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abilities behavior of the clusters with n = 19 — 22 can 
be explained in the light of their respectively geometri- 
cal characteristics. For Ge2o and Ge2i, their optimized 
structures are stacked layer configurations and have the 
comparable volume, which leads to their almost same po- 
larizabilities. In the case of Geig and Ge22 , since our cal- 
culated lowest energy structures are both near-spherical 
compact structures, they also have the tantamount po- 
larizability. 

In conclusion, we have calculated dipole moments and 
polarizabilities of Ge„(n = 2 — 25) clusters with both 
CPHF and FF approaches under B3LYP /LANL2DZ 
scheme. The main results are summarized as follow- 
ing. (1) The Gc clusters with n — 2,4,5,7 have rel- 
atively higher symmetry and the dipole moments are 
nearly zero. (2) The polarizabilities of small clusters in- 
crease rapidly in the range of 2 to 5 and fluctuate around 
5.4 A 3 /atom. Moreover, the polarizabilities of the clus- 
ters with n = 8, 11, 13, 15 are larger than the neighboring 
ones. (3) The polarizabilities are closely related to the 
HOMO-LUMO gaps and the geometrical configurations. 
The larger the HOMO-LUMO gap, the smaller the polar- 
izability of the Ge clusters. The prolate structures corre- 
sponds to relatively large polarizabilities in comparison 
with the compact structures. 
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